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Conventional computer electronics creates a di- 
chotomy between how information is processed 
and how it is stored. Silicon chips process infor- 
mation by controlling the flow of charge through 
a network of logic gates. This information is then 
stored, most commonly, by encoding it in the ori- 
entation of magnetic domains of a computer hard 
disk. The key obstacle to a more intimate integra- 
tion of magnetic materials into devices and circuit 
processing information is a lack of efficient means 
to control their magnetization. This is usually 
achieved with an external magnetic field or by the 
injection of spin-polarized currents [1, 2, 3]. The 
latter can be significantly enhanced in materials 
whose ferromagnetic properties are mediated by 
charge carriers [4]. Among these materials, con- 
ductors lacking spatial inversion symmetry cou- 
ple charge currents to spin by intrinsic spin-orbit 
(SO) interactions, inducing nonequilibrium spin 
polarization [5|, [gl, [tI, [sI, IqI, |10|, ill] tunable by local 
electric fields. Here we show that magnetization 
of a ferromagnet can be reversibly manipulated 
by the SO-induced polarization of carrier spins 
generated by unpolarized currents. Specifically, 
we demonstrate domain rotation and hysteretic 
switching of magnetization between two orthogo- 
nal easy axes in a model ferromagnetic semicon- 
ductor. 

In crystalline materials with inversion asymmetry, in- 
trinsic spin-orbit interactions (SO) couple the electron 
spin with its momentum h'k. The coupling is given 
by the Hamiltonian Tiso = t<^'^(k), where h is the 
Planck's constant and a is the electron spin operator 
(for holes a should be replaced by the total angular 
momentum J). Electron states with different sign of 
the spin projection on 11 (k) are split in energy, anal- 
ogous to the Zeeman splitting in an external magnetic 
field. In zinc-blende crystals such as GaAs there is a 
cubic Dresselhaus term[12] CId cx /c^, while strain in- 
troduces a term fl^ = CAe{kx,—ky^O) that is linear in 
/c, where As is the difference between strain in the z 
and x^y directions [13]. In wurzite crystals or in mul- 
tilayered materials with structural inversion asymmetry 
there also exists the Rashba term [14] fl^i which has a 
different symmetry with respect to the direction of /c. 



= (^R{—ky^kx,0), where z is along the axis of reduced 
symmetry. In the presence of an electric field the elec- 
trons acquire an average momentumfiAk(E), which leads 
to the generation of an electric current j = p~^E in the 
conductor, where p is the resistivity tensor. This current 
defines the preferential axis for spin precession (ri(j)). 
As a result, a nonequilibrium current-induced spin polar- 
ization (J^) II (r2(j)) is generated, whose magnitude (J^) 
depends on the strength of various mechanisms of mo- 
mentum scattering and spin relaxation[5, 15]. This spin 
polarization has been measured in non-magnetic semi- 
conductors using optical0,[l,[l,[lll,[l65 and electron spin 
resonance [17] techniques. It is convenient to parameter- 
ize (J^) in terms of an effective magnetic field H^°. Dif- 
ferent contributions to have different current depen- 
dencies ((X j or j^), as well as different symmetries with 
respect to the direction of j, as schematically shown in 
Fig. (c,d), allowing one to distinguish between spin po- 
larizations in different fields. 

In order to investigate interactions between the 
SO-generated magnetic field and magnetic domains 
we have chosen (Ga,Mn)As, a p-type ferromagnetic 
semiconductor 0, 0] with zinc blende crystalline struc- 
ture similar to GaAs. Ferromagnetic interactions in this 
material are carrier- mediated [20, 21]. The total angu- 
lar momentum of the holes J couples to the magnetic 
moment F of Mn ions via antiferromagnetic exchange 
'Hex = — • J. This interaction leads to the ferro- 
magnetic alignment of magnetic moments of Mn ions 
and equilibrium polarization of hole spins. If additional, 
non-equilibrium spin polarization of the holes (J^) is in- 
duced, the interaction of the hole spins with magnetic 
moments of Mn ions allows one to control ferromagnetism 
by manipulating J. Magnetic properties of (Ga,Mn)As 
are thus tightly related to the electronic properties of 
GaAs. For example, strain-induced spin anisotropy of 
the hole energy dispersion is largely responsible for the 
magnetic anisotropy in this material. (Ga,Mn)As, epi- 
taxially grown on (001) surface of GaAs, is compressively 
strained, which results in magnetization M lying in the 
plane of the layer perpendicular to the growth direction, 
with two easy axes along the [100] and [010] crystallo- 
graphic directions [22|, [2^. Recently, control of magneti- 
zation via strain modulation has been demonstrated [^2^. 
In this paper we use SO-generated polarization (J^) to 
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FIG. 1: Layout of the device and symmetry of the 

SO fields, a) AFM image of sample A with 8 non-magnetic 
metal contacts, b) Diagram of device orientation with re- 
spect to crystallographic axes, with easy/hard magnetization 
axes marked with blue/red dashed lines. Measured directions 
of ^ field are shown for different current directions. c,d) 
Orientation of effective magnetic field with respect to cur- 
rent direction for c) strain- induced and d) Rashba SO inter- 
actions. Current-induced Oersted field under the contacts has 
the same symmetry as the Rashba field. 



manipulate ferromagnetism. 

We report measurements on two samples fabricated 
from (Ga,Mn)As wafers with different Mn concentra- 
tions. The devices were patterned into circular islands 
with 8 non-magnetic Ohmic contacts, as shown in Fig. a 
and discussed in Methods. In the presence of a strong 
external magnetic field H, the magnetization of the fer- 
romagnetic island is aligned with the field. For weak 
fields, however, the direction of magnetization is primar- 
ily determined by magnetic anisotropy. As a small field 
(5 < < 20 mT) is rotated in the plane of the sample, 
the magnetization is re-aligned along the easy axis closest 
to the field direction. Such rotation of magnetization by 
an external field is demonstrated in Fig. . For the cur- 
rent I||[110], the measured Rxy is positive for M||[100] 
and negative for M||[010]. Note that Rxy^ and thus also 
the magnetization, switches direction when the direction 
of H is close to the hard axes [110] and [110], confirm- 
ing the cubic magnetic anisotropy of our samples. The 
switching angles (^h = ZHI where Rxy changes sign are 
denoted as (p^^ on the plot. 

In the presence of both external and SO fields, we 
expect to see a combined effect of + H on the di- 
rection of magnetization. For small currents (few /iA) 
H^^ ^ 0, and Rxy does not depend on the sign or 
the direction of the current. At large dc currents the 
value of (f^^ becomes current-dependent and we define 
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FIG. 2: Dependence of transverse anisotropic magne- 
toresistance on current and field orientation a,b) Trans- 
verse anisotropic magnetoresistance Rxy is plotted as a func- 
tion of external field direction cpn for if = 10 mT and current 
/ = ±0.7 mA in Sample A. The angles cp^^ mark magnetiza- 
tion switchings. In c) magnetization switching between [100] 
and [010] easy axes is plotted for several values of the current. 



A^^)(/) = ^^)(/) - ^^^(-/). Specifically, for I||[110] 
the switching of magnetization [010] [100] occurs for 
/ — +0.7 mA at smaller (/p^^ than for / = —0.7 mA, 
A(/p^^ < 0. For the [OlO] [100] magnetization switch- 
ing, the /-dependence of switching angle is reversed, 
Aif'^' > 0. There is no measurable difference in switch- 
ing angle for the [100] [010] and [100] [010] tran- 
sitions (A(/:?^'^^ ^ 0). When the current is rotated by 
90° (I||[110]), we observe A(/pg^ > 0, Aip^^^ < 0, and 



(1 3) (2) 

Aif fj' ^ 0. In Fig. (c) we show that A(p^^\I) decreases 
as current decreases and drops below experimental reso- 
lution of 0.5° at / < 50 /iA. Similar data is obtained for 
Sample B, see Fig. S4 in Supplementary Information. 

The data can be qualitatively understood if we con- 
sider an additional current-induced effective magnetic 
field H®^, as shown schematically in Fig. b. When an 
external field H aligns the magnetization along one of 
the hard axes, a small perpendicular field can initiate 
magnetization switching. For I||[110], the effective field 
H^^||[110] aids the [100] [010] magnetization switch- 
ing, while it hinders the [100] [010] switching. For 
(/?^^ ^ 90° and (f^ ^ 270°, where [010] ^ [100] and 
[OlO] [100] magnetization transitions occur, H®^||H 
does not affect the transition angle, A(/?^'^^ = 0. For 
I||[110] the direction of the field H^^||[110] is reversed 
relative to the direction of the current, compared to the 
I||[110] case. The symmetry of the measured H®^ with 
respect to I coincides with the unique symmetry of the 
strain-related SO field. Fig. (c). 

The dependence of Acp^^ on various magnetic fields 
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FIG. 3: Determination of current-induced effective SO 
magnetic field. a,b) Difference in switching angles for oppo- 
site current directions Acp^^ as a function of / are plotted for 
Sample A for different external fields H for orthogonal current 
directions. In c) the measured effective field H""^ = H^''±H^'' 
is plotted as a function of average current density (j) for Sam- 
ple A (triangles) and Sample B (diamonds). In d) we schemat- 
ically show different angles involved in determining H^^: (pn 
is the angle between current I and external magnetic field 
H; A(pH is the angle between total fields H + H^^(+/) and 
H + H^^(-/), and is the angle between I and H^^(+/). 



and current orientations is summarized in Fig. [3^a,b). 
Assuming that the angle of magnetization switching de- 
pends only on the total field H®^ + H, we can extract 
the magnitude H^^ and angle 6 = ZIH®^ from the mea- 
sured A(p^^ ^ thus reconstructing the whole vector H®^. 
Following a geometrical construction depicted in Fig. [3li 
and taking into account that is small, we find that 

H^^ ^ i7sin(A(^gV2)/sin(^ - (/;g^), 

and can be found from the comparison of switching at 
two angles. We find that ^ 90°, or H^^±I for /||[110] 
and /||[110]. In order to further test our procedure we 
performed similar experiments with small current / = 10 
fiA but constant additional magnetic field ^H±I playing 
the role of H®^. The measured SH{A(Ph) coincides with 
the applied SH within the precision of our measurements, 
(see Fig. S5 of Supplementary Information). 

In Fig. ^c), H^^ is plotted as a function of the av- 
erage current density {j) for both samples. There is a 
small difference in the H^^ vs (j) dependence for I||[110] 
and I||[110]. The difference can be explained by con- 
sidering the current-induced Oersted field H^^ oc / in 
the metal contacts. The Oersted field is localized under 
the pads, which constitutes only 7% (2.5%) of the total 
area for samples A (B). The Oersted field has the sym- 
metry of the field depicted in Fig. (d), and is added to 
or subtracted from the SO field, depending on the cur- 
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FIG. 4: Current-induced reversible magnetization 
switching a) dependence of Rxy near the [010] [100] 
magnetization switching is plotted for / = ±0.7 mA in Sample 
A for /||[110]. h)Rxy shows hysteresis as a function of current 
for a fixed field H — Q mT applied at LpH — 72°. c) Mag- 
netization switches between [010] and [100] directions when 
alternating ±1.0 mA current pulses are applied. The pulses 
have 100ms duration and are shown schematically above the 
data curve. Rxy is measured with / = 10 /iA. 



rent direction. Thus, H^^ = + H^^ for I||[110] and 
^eff ^ H^o_jjOe j|| [;l];q] estimate the fields to be 

as high as 0.6 mT under the contacts at / = 1 mA, which 
corresponds to H^^ ^ 0.04 mT (0.015 mT) averaged over 
the sample area for samples A (B). These estimates are 
reasonably consistent with the measured values of 0.07 
mT (0.03 mT). Finally, we determine H^^ as an average 
of H^^ between the two current directions. The SO field 
depends linearly on j, as expected for strain-related SO 
interactions: dH'^'/dj = 0.53 • 10"^ and 0.23 • 10"^ T 
cm^/A for samples A and B respectively. 

We now compare the experimentally measured 
with theoretically calculated effective SO field. In 
(Ga,Mn)As, the only term allowed by symmetry that 
generates H^^ linear in the electric current is the term, 
which results in the directional dependence of on j 
precisely as observed in experiment. As for the mag- 
nitude of for three-dimensional J = 3/2 holes we 
obtain 



H^°(E) = 



eCAs {-SSuhTh + ISmn) 



{Ex, —Ey, 0), 



where E is the electric field, g* is the Luttinger Lande 
factor for holes, /i^ is the Bohr magneton, and rih^i and 
Th^i are densities and lifetimes for the heavy (h) and light 
(1) holes. Detailed derivation of H^^ is given in Sup- 
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plementary Information. Using this result, we estimate 
dH'^'/dj = 0.6 • 10"^ T cmVA assuming rih = n -> ni 
and Th = mh/{e^pn)^ where p is the resistivity mea- 
sured experimentahy, and using = 10~^, n = 2 • 10^^ 
cm~^. The agreement between theory and experiment 
is excehent. It is important to note, though, that we 
used GaAs band parameters [25| rrih = 0.4mo, where mo 
is the free electron mass, 0* = 1.2 and C = 2.1 eV-A. 
While the corresponding parameters for (Ga,Mn)As are 
not known, the use of GaAs parameters appears reason- 
able. We note, for example, that GaAs parameters ad- 
equately described tunnelling anisotropic magnetoresis- 
tance in recent experiments [26]. 

Finally, we demonstrate that the current-induced ef- 
fective SO field H^^ is sufficient to reversibly manipu- 
late the direction of magnetization. In Fig. [4^ we plot 
the (/:?i:/-dependence of Rxy for Sample A, showing the 
[010] [100] magnetization switching. If we fix = 6 
mT at cpH = 72°, R^y forms a hysteresis loop as cur- 
rent is swept between ±1 mA. Rxy is changing between 
±5 1^, indicating that M is switching between [010] and 
[100] directions. Short (100 msec) 1 mA current pulses of 
alternating polarity are sufficient to permanently rotate 
the direction of magnetization. The device thus performs 
as a non-volatile memory cell, with two states encoded 
in the magnetization direction, the direction being con- 
trolled by the unpolarized current passing through the 
device. The device can be potentially operated as a 4- 
state memory cell if both [110] and [110] directions can 
be used to inject current. We ffnd that we can reversibly 
switch the magnetization with currents as low as 0.5 mA 
(current densities 7 • 10^ A/cm^), an order of magnitude 
smaller than by polarized current injection in ferromag- 
netic metals [l|, H, [sj, and just a few times larger than 
by externally polarized current injection in ferromagnetic 
semiconductors Q . 

Methods 

The (Ga,Mn)As wafers were grown by molecular beam 
epitaxy at 265 °C and subsequently annealed at 280 
°C for 1 hour in nitrogen atmosphere. Sample A was 
fabricated from 15-nm thick epilayer with 6%Mn, and 
Sample B from 10-nm epilayer with 7% Mn. Both 
wafers have Curie temperature Tc ~ 80 K. The de- 
vices were patterned into 6 and 10 /im-diameter cir- 
cular islands in order to decrease domain pinning. 
Cr/Zn/Au (5nm/10nm/300nm) Ohmic contacts were 
thermally evaporated. All measurements were performed 
in a variable temperature cryostat at T = 40 K for Sam- 
ple A and at 25 K for Sample B, well below the temper- 
ature of (Ga,Mn)As-speciffc cubic-to-uniaxial magnetic 
anisotropy transitions [27], which has been measured to 
be at 60 K and 50 K for the two wafers. Temperature 
rise for the largest currents used in the reported experi- 
ments was measured to be < 3 K. 

Transverse anisotropic magnetoresistance Rxy = 
Vy/Ix is measured using the four-probe technique, which 



insures that possible interfacial resistances, e.g., those 
related to the antiferromagnetic ordering in the Cr wet- 
ting layer [i^, do not contribute to the measured Rxy 
The DC current Ix was applied either along [110] (con- 
tacts 4-8 in Fig. a) or along [lIO] (contacts 2-6) direction. 
Transverse voltage was measured in the Hall conffgura- 
tion, e.g., between contacts 2-6 for /a^||[110]. To ensure 
uniform magnetization of the island, magnetic ffeld was 
ramped to 0.5 T after adjusting of the current at the 
beginning of each ffeld rotation scan. We monitor Vx 
between different contact sets (e.g. 1-7, 4-6 and 3-5) to 
conffrm the uniformity of magnetization within the is- 
land. 

In order to determine the direction of magnetization 
M, we use the dependence of Rxy on magnetization [ioj : 

Rxy Ap sin (fM cos (fM , 

where Ap = p\\—p±^ p\\ < p± are the resistivities for mag- 
netization oriented parallel and perpendicular to the cur- 
rent, and (fM = ZMI is an angle between magnetization 
and current. In a circular sample the current distribu- 
tion is non-uniform and the angle between the magneti- 
zation and the local current density varies throughout the 
sample. However, the resulting transverse AMR depends 
only on cpM- For the current-to-current-density conver- 
sion, we model our sample as a perfect disc with two point 
contacts across the diameter. The average current den- 
sity in the direction of current injection is (j) = 21/ {nad)^ 
where a is the disk radius and d is the (Ga,Mn)As layer 
thickness. In a real sample the length of contact over- 
lap with (Ga,Mn)As insures that j changes by less than 
factor of 3 throughout the sample. A detailed discussion 
of the current distribution and of measurements of Joule 
heating can be found in Supplementary Information. 
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JOULE HEATING 



(Ga,Mn)As is a magnetic semiconductor with strong temperature dependence of resistivity, see Fig. [51(a). The 
enhancement of resistivity at 80 K is due to the enhancement of spin scattering in the vicinity of the Curie temperature 
Tc. Inelastic scattering length in these materials is just a few tenths of nm, and we expect holes to be in thermal 
equilibrium with the lattice[30[. Thus resistivity can be used to measure the temperature of the sample. 




FIG. 5: Current-induced heating a) Temperature dependence of resistivity for sample A; b) current and c) temperature 
dependence of sample resistance in the vicinity of 40 K; d) sample heating as a function of dc current. 



In Fig. [5fb,c) we plot temperature and current dependences of the sample resistance in the vicinity of 40 K. This 
data is combined in (d), where the sample temperature change AT due to Joule heating is plotted as a function of dc 
current. The maximum temperature rise does not exceed 3 K at / = 0.7 mA in our experiments. This small heating 
ensures that the sample temperature stays well below the Curie temperature (~ 80 K) and the (Ga,Mn)-specific 
cubic-to-uniaxial magnetic anisotropy transition 60 K for sample A and ^ 50 K for sample B) when experiments 
are performed at 40 K and 25 K for samples A and B, respectively. Observation of different angles for magnetization 
switching for +/ and — / (Fig. 2) further confirms that heating is not responsible for the reported effects (Joule heating 
is cx and does not depend on the current direction). 



CURRENT DISTRIBUTION IN CIRCULAR SAMPLES 



Magnetization-dependent scattering in (Ga,Mn)As results in an anisotropic correction to the resistivity tensor p 
which depends on the angle cpm between magnetization M and local current density j [29,] : 



Pxx = P± + {P\\ - P±) COS^((/Prn), 
Pxy = {P\\ - P±) Sm{(pm) COs{(pm), 



(1) 
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where p\\ {p±) are the resistivities for j||M (j±M), and we assumed that both j and M he within the plane of the 
sample. The off-diagonal resistivity (transverse anisotropic magnetoresistance) pxy can be non-zero even in the absence 
of the external magnetic field. The difference {p\\ — p±)/p± ^ 0.01 and we first calculate the local potential <po{x^y) 
inside the sample by approximating it as a disk of radius a and thickness d with isotropic resistivity po = (py + p^)/2: 



ird 



In 



(a - x)^ + y'^ 



{a + + y'^ 



(2) 



where current / is injected along the x-axis. Current density j = V0o/po is plotted in Fig. [BJa). Metal contacts have 
a radius of ~ 0.5 pm in our samples, which limits the current density near the current injection regions. Integrating 
j over the sample area we find average current density 



(3) 




FIG. 6: Current distribution a) Vector plot of local current density distribution in the sample; b) angles between 

magnetization M and current I||x are defined; c) Color map plot of Oersted field (i^;^) distribution in a disk-shaped 

sample. 

We find the transverse voltage Vy as a correction to the (j)o potential due to the anisotropic resistivity p\\ — p± 0: 

/ao 
[-pxy- jx{y) + pxx • jy{y)]dy. (4) 
-ao 

The current distribution is non-uniform, and the local electric field depends on the total angle cpm = — , where 
(fM = MI and cpj = jl, see Fig. [6l^b). This integral can be evaluated analytically, and the transverse anisotropic 
magnetoresistance (AMR) Rxy is found to be the same as for an isotropic current fiow, independent of the distance 
xo of the voltage contacts from the center of the disk: 

Rxy = Vy/I = {p\\ - pi_) cos{(Pm) sm{(pM)' (5) 
The magnetization angle cpM can therefore be directly calculated from the measured transverse resistance Rxy. 



CURRENT-GENERATED OERSTED MAGNETIC FIELDS 



In this section we estimate conventional current-generate magnetic fields in our device that are not related to spin- 
orbit interactions. There are two contributions to the Oersted magnetic fields: a magnetic field due to non-uniform 
current distribution within the sample, and a field generated by high currents in the vicinity of the metal contacts. 

We can calculate the Oersted field inside (Ga,Mn)As by using the Biot-Savart formula: 



(6) 
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spin-orbit field contact Oersted field 

FIG. 7: Oersted field a) schematic illustration of the origin of the in-plane Oersted field H^^ under gold contact pads; b,c) 
symmetry of H^^ and H^^ fields. 

where /io is the permeability of free space, and the integral is taken over the volume of the disk. The most significant 
H-^^ normal component of the field is shown in Fig.[3l^c). The largest H-^^ < 1 Oe, which is negligible compared to the 
2000 Oe anisotropy field that keeps the magnetization in-plane. 

The second contribution to the Oersted field originates from contact pads, see Fig. [71 The conductivity of gold 
contacts is much higher that of (Ga,Mn)As, and the current flows predominantly through the metal within contact 
regions, thus generating both in-plane {H^^) and out-of plane (H^^) magnetic fields in (Ga,Mn)As underneath and 

at the edges of the contact pads. The maximum value of the field can be estimated as H^^ ^ H^^ = fioI/2w, where 
/ is the total current and w = 1 /im is the width of the contact pad. This field can be as high as 6 Oe for / = 1 mA. 
The field is localized under the pads, which constitute only l/12th of the sample area. 

The H^^ field does not induce in-plane magnetization rotation. The H^^ field and the effective spin-orbit field have 
different symmetries with respect to the current rotation, see Fig.[3fb,c), and thus can be experimentally distinguished. 
The two fields point in the same direction for I||[110], but in the opposite direction for the current rotated by 90°, 
I||[ll0]. Experimentally, we observe an effective field which corresponds to the symmetry of the SO effective field. 
However, there is a small difference in the slopes of A(j)H vs / curves for the two orthogonal current directions. 
Fig. 3(a,b), because the contact field is added to the SO field for I||[110] and subtracted from SO field for I||[110]. 
Both fields ex /. From the ratio of the slopes (~ 1.2) we can calculate the strength of the contact field, H^^ ^ O.IH^^. 
This experimentally found ratio is consistent with the above estimate if we average the contact Oersted field over the 
sample area. 



DEPENDENCE OF TRANSVERSE ANISOTROPIC MAGNETORESISTANCE ON CURRENT AND 

FIELD ORIENTATION 




FIG. 8: Dependence of transverse anisotropic magnetoresistance on current and field orientation for Sample 

B. Transverse anisotropic magnetoresistance Rxy is plotted for current I||[110] (a) and I||[110] (b) for / = ±0.75 mA with 
constant magnetic field H — 20 mT as a function of field angle (fH- 



In order to test the procedure of the effective field mapping, we performed control experiments where a small 
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constant external magnetic field ^H±I was playing the role of spin-orbit field. In these experiments the current was 
reduced to / = 10 /iA. The results, shown in Fig. [9l are quantitatively similar to the effect of the spin-orbit field. 




FIG. 9: Control experiment with additional external field a) Rxy is plotted for Sample A with Htotai = ^H + H, where 
SH1.I and / = 10 /iA, 6H = 0, 0.8 mT. For comparison, in b) similar data are plotted for SH = but / = ±0.5 mA. 



CALCULATION OF THE EFFECTIVE SPIN-ORBIT FIELD INDUCED BY THE ELECTRIC CURRENT 

Manipulation of localized spins, electronic, nuclear or ionic, can be achieved via manipulation of free carrier spins. 
The free carrier spins can be manipulated by the external magnetic field, by the Oersted magnetic field of the 
current, and by the electric current via intrinsic spin-orbit interactions. The intrinsic spin-orbit interactions arise in 
crystalline systems, in which axial vectors, such as spin polarization, and polar vectors, such as the electric current, 
behave equivalently with respect to the symmetry trasformations of a crystal. The crystal symmetry then allows 
the transformation of the electric current into a spin polarization of charge carriers. In this work, Mn ions of the 
ferromagnetic semiconductor (Ga,Mn)As, and thus its ferromagnetic properties, are affected by the electric current 
via the intrinsic spin-orbit interactions. 

In (Ga,Mn)As, charge carriers are holes with an angular momentum J = 3/2. In contrast to electron systems, the 
hole system is defined by a very strong coupling of the total angular momentum J to the hole momentum p, which 
includes both terms quadratic in p and independent of p. These terms are quadratic in J, and they are not present 
for electrons with spin 1/2. The Luttinger-Pikus Hamiltonian quadratic in J is[13] 

Hh = Aop^ + ^1 XI '^iP^^ + ^2 JiJjPiPj + ^1 XI ^^^^^ + ^2 X ^^^i^^i' C^) 

where i, j = x, z. Despite the presence of a very strong spin-orbit coupling, which leads to the spectral splitting of 
holes into two pairs of states, this Hamiltonian on its own cannot result in a spin polarization of holes induced by 
the electric current. Terms capable of generating spin polarization in systems characterized by the absence of center 
of symmetry in the crystal and by a corresponding additional lowering of the crystalline symmetry in the presence of 
strain, read 

= iv^JiPi{pl+i - pl+2) ^ c^[JiPi{£iJ,i^iJ,i -£i+2,i+2) + {JiPi+1 - (8) 

i i 

where cyclic permutation of indices is implied. The first term is cubic in the hole momentum, and it can lead only 
to the polarization of hole spins cubic in the electric current (and only when the current direction is away from the 
high symmetry axes). For effects linear in electric current this term is only relevant insofar as it contributes to the 
spin relaxation of the holes. The third term contains off-diagonal components of the strain tensor, and is negligible 
in (Ga,Mn)As crystals under study. In this system, strain originates from doping by Mn ions, and constitutes 
tension along the growth axis 2:||[001] defined by the component Szz and As = Szz — ^xx = ^zz — ^yy Thus only 
the second term results in a current-induced spin polarization. The symmetry of the corresponding effective field, 
r2(p) = CAe{px, ~Py^ 0)5 depends markedly on the crystallographic orientation. When an electric field is applied, the 
direction of the generated hole spin polarization with respect to the orientation of the electric current is the same as 



10 



the direction of the SO effective field with respect to the hole momentum. Such peculiar symmetry differs from the 
symmetry of the Oersted magnetic field, and thus allows one to distinguish between these effects. 

We consider now the approximation linear in strain, when only the strain-dependent term proportional to C is 
taken into account, and strain-dependent terms in TL^ are omitted. In this case the hole spectrum given by Tih splits 
into heavy (h) and light (/) hole branches. The mechanism of generation a spin polarization by the effective SO 
field in the presence of an electric current is simply a shift in the distribution functions for heavy and light holes in 
momentum space. In contrast low symmetry electron systems |15|] , where spin polarization is associated entirely with 
the relaxation of spins, in case of holes the spin relaxation occurs on the time scale of momentum relaxation and plays 
no role in the current-induced spin polarization. At low temperatures the hole angular momentum density is given by 



where z = 1, 2 correspond to principal axes x and characteristic times Th^i are defined by mobilities of holes in the 
corresponding bands, and ny^^^^ are densities of holes in these bands. At room temperatures Ep in the denominator 
is to be replaced by 3/2/c^T, T being the lattice temperature and the Boltzman constant. Estimates show that 
the negative term in brackets of Eq. [9]is dominant. 

We note that in the case of very strong deformations the spin relaxation of holes occurs on the times scale longer 
than that of momentum relaxation. Then simple shift of hole distribution functions in momentum space is no longer 
sufficient for generating spin polarization by current, and the mechanism of the effect becomes analogous to that for 
electrons [15]. We will present the results for hole spin polarization generated by electric current at arbitrary value of 
strain elsewhere. 

The spin polarization given by Eq. [9] leads to an effective magnetic field acting on the Mn ions. In order to calculate 
what external magnetic field would result in the same polarization as that generated by the current, we calculate the 
average spin density induced by an external magnetic field: 



The ratio of polarizations {J^ ^) and (J^^ ^) gives the electric field polarization measured in units of magnetic field. 
We note that while the SO field affects Mn ions only via the exchange interaction, the Oersted or the external magnetic 
field also acts on the ions directly. However, the magnitude of the exchange interaction, A = — 5 meV is quite large, 
making the exchange interaction dominant. We will therefore omit the discussion of direct polarization of Mn by 
external fields. 




(9) 



■(H) V _ ilg*HBH{nh + ni) 



(10) 



